Infinite Invariant Density Determines Statistics of Time Averages for Weak Chaos 
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Weakly chaotic non-linear maps with marginal fixed points have an infinite invariant measure. 
Time averages of integrable and non-integrable observables remain random even in the long time 
limit. Temporal averages of integrable observables are described by the Aaronson-Darling-Kac 
theorem. We find the distribution of time averages of non-integrable observables, for example the 
time average position of the particle, x. We show how this distribution is related to the infinite 
invariant density. We establish four identities between amplitude ratios controlling the statistics of 
the problem. 
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Low dimensional chaotic systems by definition have 
positive Lyapunov exponents and have been extensively 
used to test basic assumptions of statistical physics. 
Weakly chaotic systems have zero Lyapunov exponents, 
namely the separation of trajectories is sub-exponential, 
though the deterministic motion remains quasi-random. 
In many cases discrete maps are used to model the dy- 
namics, since they help to establish a deep understanding 
of the fundamental issues without being too complicated 
(importantly numerics converge faster than in more real- 
istic models). In particular Pomeau-Manneville [if maps 
are weakly chaotic 0] and are characterized by marginal 
instability. These maps were used to model intermittency 
anomalous diffusion and aging [?[. Such sys- 

tems are described by an infinite invariant density (ooD) 
[1, Q: a non-normalizable density defined below. It is 
also well known that temporal averages in such systems 
are not equal to a corresponding ensemble average, in- 
stead time averages remain random variables even in the 
long measurement time limit [i[iMl. Since chaos is a 
precondition for statistical physics, it is not very surpris- 
ing that weak chaos implies the breakdown of standard 
ergodic theory. 

For an ergodic process, in the long time limit the time 
average of an observable is equal to the corresponding 
ensemble average. The ensemble average and hence the 
time average are obtained from the normalized invari- 
ant density, if it exists. A fundamental extension of 
standard ergodic theory is to find the distribution of 
time averages of generic observables for weakly chaotic 
systems where the underlying invariant density is non- 
normalizable. The Aaronson-Darling-Kac (ADK) theo- 
rem Q gives a partial answer to this problem. Briefly, 
an observable whose average with respect to the ooD is 
finite, the distribution of properly scaled temporal aver- 
ages is the Mittag-Lefher distribution. The ooD is essen- 
tial for the description of these fluctuations. For example, 
the separation of trajectories is described by a stretched 
exponential (a manifestation of weak chaos) and the dis- 
tribution of separation rates is provided by the ADK the- 
orem [illil. 

In this manuscript, we consider the very large class of 
non-integrable observables. We focus on the position of 




FIG. 1: (Color online) (a) The map Eq. ^ with L = 3 
and a — 3/4 has IFPs on 0, 1, 2. (b) The ooD exhibits non- 
integrable divergence on the IFPs (parameters are given in 
third example in the text). 



a particle xt in an interval (0,i) and obtain the distri- 
bution of its time average. Importantly we show how 
the distribution of time averages of non integrable ob- 
servables is related to the underlying ooD. Previously 
Thaler and Zweimiiller [lol [Tl| considered an important 
non-integrable observable: the occupation fraction; i.e., 
the fraction of time the particle spends within a given 
domain. They rigorously showed it is described by the 
Lamperti distribution (see details below). We provide a 
very general conjecture for the distribution of time av- 
erages of non-integrable observables, without giving a 
rigorous proof but rather relying on simple arguments. 
Further we derive the identity of four amplitude ratios 
which govern the statistics of the problem. These identi- 
ties bridge between the stochastic and dynamical theories 
in this field. 

Model and observable. We consider measure preserv- 
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ing maps Xt+i = M{xt) with xt € (0,-L). Our observ- 
able is Xt and our goal is to calculate the distribution 
of its time average x = X]t=o^t/^i the limit of long 
time. We assume that the map has N indifferent fixed 
points (IFPs) located on {a;(l), ...x{j), ...x{Ny\ such that 
M{x) ~ x + 2^/°'aj\x-x{j)\^+^/°' asx^ x{j) and ^ 
(IFPs are also called marginal fixed points). Throughout 
this work, j is a label of the IFPs. Here we consider 
< a < 1 since in that regime the distribution of x is 
non-trivial. An example map is shown in Fig. [TJ We con- 
sider maps where the trajectory of the particle visits the 
vicinity of all the IFPs; i.e., we exclude stable points or 
a decomposable phase space, so the transformation has 
an infinite invariant measure. Such maps exhibit non- 
Gaussian intermittency and hence have attracted vast 
research using various methods such as CTRW and 
periodic orbit theory Q. 

Power law sojourn times are related to the injection 
probability. Let us consider the IFP a;(l) which we des- 
ignate to be on the origin x(l) = 0. In the vicinity of 
this point the map is xt+i ~ -I- 2^/°'ai{xt)^/°'~^^ for 

> and < a < 1 while ai > (other IFPs have con- 
stants Oj). Starting on xq the time r it takes the particle 
to reach a threshold Xc is determined by the continuous 
approximation of the map dx/dt ~ 2^/"ai(a::)^/"+^. Fol- 
lowing Geisel and co-workers P| this gives 
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During the evolution the particle is injected in the vicin- 
ity of the IFP many times and hence xq is treated 
as a random variable whose probability density func- 
tion (PDF) is P™[xo]. It follows that the waiting time 
r, the time the particle remains in the vicinity of the 
j ~ 1 IFP, is a random variable with the PDF iPi{t) — 
P"^[xo]\dxo/dT\. A similar formula holds for the jth IFP. 
Using Eq. ([T]) one finds the PDF of waiting times Q 



V'j(t) - y4jT-(i+") with Aj = P'" [x{j)] 
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Notice that this result is independent of the threshold 
Xc- Here it is assumed that the injection PDF P™[a;(j)] 
is smooth in the vicinity of the IFP. Eq. ([2]) is well 
known but actually rather formal since it expresses ipj (r) 
in terms of the unknown injection PDF. Below we will 
relate the injection PDF with the ooD. The power-law 
PDF Eq. ([2]) indicates a diverging average sojourn time 
since < a < 1. The corresponding stochastic picture 
[1, § is of a particle jumping between neighborhoods of 
the IFPs {x{l) ■ ■ ■ x{N)} with power law sojourn times 
for the trapping events. 

The infinite invariant density plays a crucial role and 
it is defined as flql 



p{x)~p{x,t)/t''- 



(3) 



where p{x, t) is the density of particles [in simulations 
we use initial conditions uniformly spread in (0,L)]. 



When a < 1, the invariant density is non-normalizable, 
'p{x)dx — oo, and hence its name. Such ooDs are 
not common in physics though recently an application 
was suggested in the context of cooling in optical lattices 
[13] ■ Note that the density p{x,t) is, as usual, normahz- 
able for any finite t since the maps conserve the number 
of particles. In the vicinity of the IFP x{j) one finds the 
non-integrable behavior 



p{x) 



-I/O 
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where bj > is an amplitude which is generally unknown. 
An example ooD is shown in Fig. [1] based on a numerical 
simulation which allows us to estimate the bjS. 

To understand better such a behavior we use simple 
arguments. Note that the density normalized to unity is 



p{x, t)dx ~ t 



It 



(5) 



where tx^x+Ax is the time the particle spends in (a:, x^dx) 
p^ . Let us consider the vicinity of the first IFP x(l) = 0. 
The time tx,x+Ax is proportional to Nr: the number 
of injections to the vicinity of the IFPs, multiplied by 
P'"[a;(l)]da; [which gives the number of visits in the in- 
terval (a;, x~\-dxy\. t^^x+dx is also proportional to the time 
the particle stays in (x, x + dx) during each visit, which 
we call M. Thus, close to the IFP, 



p{x,t) 



iVflP'" [x(l)] At 



(6) 



As is well known from renewal theory Q the num- 
ber of renewals or injections scales like Nr ~ Ct". The 
pre-factor C can be roughly estimated however below we 
show that it is an irrelevant parameter. Using Eq. ([T]) 
we have when x — > 



At: 
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so that the closer the particle is to the IFP x(l) = the 
longer is At. Similarly we analyze other IFPs. Putting 
these pieces of information together, we find 



p{x,t) 



4^where5, = ^^^^;iM)l, 



(8) 

Eq. ([8]) shows, a relation between the amplitudes of the 
ooDs; i.e., the bjS and the injection probabilities P™[x(j)] 

The time average x is now considered. During the 
evolution the trajectory of the particle Xt spends long 
times, of the order of the measurement time in the vicin- 
ity of the IFPs. In contrast the time it takes the particle 
to jump between one IFP state to another is short and 
can be neglected. Hence along a trajectory xt attains ob- 
servable values which are (nearly) equal to the locations 
of the IFP {a;(l) • • ■x{N)}. In each one of these states 
the particle remains a time tj with j = 1 , • ■ ■ ^ N which 
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FIG. 2: (Color online) Numerical PDF of ^ (solid line), per- 
fectly matches the theory Eq. (|20|l (dashed line) without any 
fitting [map Eq. (|19[) . left panel a = 0.3, right a — 0.75 and 



is the occupation time of state j. It follows that the time 
average is 



t 



(9) 



Each tj is a sum of many independent identically dis- 
tributed random sojourn times drawn from the long 
tailed PDF i/jjir), Eq. ([2|). Hence the occupation time tj 
is distributed according to Levy statistics, i.e. the gen- 
eralized central limit theorem holds. More precisely tj 
is a stable random variable whose PDF is the one sided 
Levy function with index < a < 1. Let p^'^ — {tj)/t 
be the averaged occupation fraction treated rigorously in 
[Tol [TH which is nothing but the probability that a mem- 
ber of an ensemble of noninteracting particles is in the 
vicinity of the IFP j. Since the occupation time tj scale 
with Aj and t — S^* 



eq 



A, 



(10) 



where Aj is the amplitude of the waiting time PDF, de- 
fined in Eq. Importantly, using Eqs. (|2ll0p . 



[a;(j)]|a,|- 



(11) 



which relates occupation fractions with injection proba- 
bilities. Using Eqs. ((STTTt 



cq _ "Jl^-Jl 
Pj ~ 



-a+1 



Ei=ibj\aj\ 



(12) 



which relates the occupation fractions and the ooD. 

The distribution of observables like x was recently 
studied within the continuous time random walk model, 
a stochastic approach extensively applied, though so far 
without an underlying ooD. Briefly, as mentioned, tj is a 
stable random variable, and since x [Eq. ^] is a linear 



combination of such independent random variables, one 
finds the PDF of the time average (20j 



fa (x) = lim Im^ 



cq I — 



(13) 



where i = v^— 1 and Im denotes the imaginary part. We 
see that the PDF of x is controlled by the nonlinear- 
ity of the map in the vicinity of the IFPs, i.e. a, the 
values of the observable on these points {a;(j)}, the equi- 
librium probabilities p^** which in turn depend on either 
the ooD, Eq. or the injection PDF, Eq. Thus 

once the invariant density is known one may obtain full 
information on the fluctuations of the time average of 
our observable. Notice that when a — > 1, Eq. (fTS]) 
yields limc^i [x) ~ 6{x - (x)) where {x) = J2pT^^^^ 
is the ensemble average. For a general non-integrable 
observable 0{xt), the distribution of the time average 
O = J2lzl 0(xt)/t is fa (O) as in Eq. ^ where on the 
right hand side we replace x{j) with 0{x{j)). 

A first illustration will be a system with two IFPs. We 
consider xt £ (0, 1) and 

r a:t + 2i/"(a;t)i+i/" 0<xt<l/5 
M{xt) = < 1 + ^rm^ 1/5 < < 11/20 

[ Xt-2^^°'{l-Xt)^+^/°' 11/20 <2:t<l, 

(14) 

hence x(l) = and x(2) = 1 are the IFPs of the map 
and |ai| — [02 1. We first concentrate on the injection 
PDF P'"[a;]. We partition the map into two parts with 
a boundary on < < 1. Following a trajectory we 
record events where the particle jumps over the bound- 
ary, either from left to right or vice versa. Each time 
the particle is injected into one of the domains a; < a;c or 
X > a;c we record its landing position and thus generate 
a histogram which gives P™[x]. Not surprisingly P'"[x] 
will depend on the choice of Xc- However, interestingly, 
the ratio P™[a:;(l)]/P™[a;(2)] is a constant independent of 
the value of Xc- To understand this behavior note that 
according to Eq. ([8]) we get the amplitude ratio relation 



= P'"' [xi2)] /P- [:e(1)] 



(15) 



and since 62/61 is clearly Xc independent so is the right 
hand side of this equation. Starting with a uniform den- 
sity we evolve the system until time 10"*, obtain an es- 
timate for the ooD p(a::), and with it find 61 and 62. 
For a = 0.75 we find 61 = 0.075,62 = 0.16 and for 
= 0.5 P'"[x(l)] = 0.86 and P'"[a;(2)] = 1.86 while 
P"^[x{l)] = 1.18 and P'"[2;(2)] = 2.58 for x^ = 0.3. 
Hence Eq. p5|) stands the numerical test. We have also 
verified this equation with other parameters. 

After we get the amplitudes of the infinite invariant 
density, 61 and 62, we may calculate pl^ and p'^'^ and so 
using Eq. ([13]) we find the PDF of x 



fa (x) 



TT ^ sin (tto) TZx°' ^ {1 — x) 



a-l 



7^2 (1 - xf" + ixf" + 27^ cos nail- xf x^ ' 

(16) 
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FIG. 3: (Color online) Numerical simulations give the PDF 
of X (solid line) that matches the analytical density (dashed 
line) Eq. ^ [map Eq. (US]), t = 10*^, L = 3 and q = 3/4]. 



which is the Lamperti PDF. The same distribution was 
previously obtained for the occupation fraction [lol - [l2j . 
As pointed out by Akimoto [13] this is not surprising 
since both observables are identical on the IFPs [for the 
occupation fraction the observable is the step function 
which is 1 on a; (2) — 1 and zero on a;(l) = 0]. The 
parameter TZ is 



^ _ ^2 _ P'"[x(2)] 



bi 



cq 
P2 



= ?-^- (17) 



cq 
Pi 



Hence one has four amplitude ratios related to the wait- 
ing times, the injection probabilities, the ooD and the 
population probabilities which determine the PDF of x. 
Amplitude ratios can be easily generalized for the case of 
iV IFPs and for the case where the OjS are not all equal 



A, 



|a,|-"F-[x(j)] _ b,\a,\- 



|a,|-"Pi"[x(z)] 6,|a,|-"+i 



(18) 



The second illustration concerns maps with N degen- 
erate IFPs. We consider N — 2L with L ~ 8 and 
xt G (-1/2,7.5). The map is 



M{xt) = xt+ 

where fc = 0, 
0-,a;(2) = 0^ 



2^/"ak{xt - fc)i+i/" 



k < Xt < k + 1/2 
k - 1/2 < Xt < k, 
(19) 

,L-1. Here 16 IFPs are on {a;(l) = 
..a;(15) = 7-,x{16) = 7+}. We use pe- 
riodic boundary conditions: if xt > 7.5 or xt < —1/2 
we transform xt to Xt — 8 ot xt + 8 respectively. We 
set all flfc = 1. Then from symmetry we expect that all 



the amplitudes bj will be identical. It then follows that 
P'i'^ = 1/{2L). For this degenerate case we get 



fa (S) = lim Im V — . a ■ 

j=o ^ -J ' ' 



(20) 



Thus due to symmetry the distribution of x depends only 
on a single parameter which is a. In Fig. [5] we show 
the PDF of X obtained numerically together with theory 
Eq. ((20)) . For a = 0.3 the distribution is wider than 
the case a = 3/4 since we expect that as q ^ 1 the 
fluctuations will vanish. Notice that fa{x) diverges on 
the IFPs reflecting trajectories with a trapping time of 
the order of the measurement time on one of these points. 

The third example is the map Eq. (|19p with L — 3 
and hence IFPs are on a:(l) = 0~,x(2) = 0"'",x(3) = 
l-,a;(4) = l+,a;(5) = 2", a;(6) = 2+. It follows that in 
the long time limit x S (0, 2). We use di = 1.1, 0,2 — 1.5 
and 03 = 2.1. We numerically obtain the ooD for a = 3/4 
and estimate 61 = 62 — 0.059, 63 64 ~ 0.04, 65 = feg — 
0.018. Inserting these values in Eq. (IT2|) we find: p^"^ = 
pI"" = 0.239, pI"" = pI^ = 0.175 and p"^ ^ pi'' = 0.086. 
This is compared with direct numerical computation of 
the occupation fraction: p'^^ — p'i^ ~ 0.242, pg'^ — p'^ ~ 
0.167,^5'' — ~ 0.091. Deviations between the two 
methods are related to the divergence of the ooD next to 
IFPs which induces errors in the estimation of the bjS. 
Inserting the latter values of pj's into Eq. we obtain 
the PDF of X which as shown in Fig. [3] perfectly matches 
direct numerical simulation. 

Discussion. We obtained the distribution of time aver- 
ages of non-integrable observables for systems with IFPs 
with an infinite invariant measure. The ooD^ the occu- 
pation fractions, the injection probabilities, and the am- 
plitudes Aj of the scale free distributions of the sojourn 
times, are all related and can be used to determine the 
non-trivial distribution of the temporal averages. There 
exists a vast number of physical systems with dynamics 
governed by power law trapping times similar to the maps 
under investigation. A fundamental experimental ques- 
tion is whether such systems, e.g. blinkirig quantum dots 
PH , two dimensional rotating flows [2I, [2^ and electro- 
hydrodynamic convection in liquid crystals [IJ] possess 
an infinite invariant measure. Hence it would be inter- 
esting to extract the invariant density from the trajecto- 
ries in such experiments. If it is of infinite measure, one 
could then use our theory to predict the distribution of 
the temporal averages. 
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